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Abstract. 

The large supercluster structures obvious in recent galaxy redshift 

■ surveys are quantified using an one-dimensional cluster analysis (core 

sampling) and a three-dimensional cluster analysis based on the minimal 
spanning tree. The comparison with the LCRS reveals promising stable 
results. At a mean overdensity of about ten, the supercluster systems 
form huge wall- like structures comprising about 40% of all galaxies. The 
overdense clusters have a low mean transverse velocity dispersion of about 
400 km/s, i.e. they look quite narrow in redshift space. We performed 
N-body simulations with large box sizes for six cosmological scenarios. 

j The quantitative analysis shows that the observed structures can be un- 

derstood best in low density models with Q m < 0.5 with or without a 
cosmological constant. 



a: 
6 

> 



1. Introduction 

Very large structures in the universe as the great attractor (Dressier et al. 1988) 
or the great wall (Ramella, Geller, & Huchra 1992) were firstly considered as 
rare peculiarities in the large scale galaxy distribution. But recently published 
rich galaxy redshift catalogues as the Durham/UKST Galaxy Redshift Survey 
(Ratcliffe et al. 1996) and the Las Campanas Redshift Survey (LCRS) (Schect- 
man et al. 1996) have established the existence of similar large structures as a 
typical characteristic of the galaxy distribution. They correspond to the super- 
clusters first identified in the distribution of clusters of galaxies (Abell, 1958; 
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Figure 1. Wedge diagram of 80° width and 1.5° depth of a mock 
catalog for the magnitude limited LCRS (7500 galaxies with 15 < m < 
18.) 

Oort, 1963), and we shall call them in the following simply superclusters. These 
structures exhibit a clear filamentary network, and the filaments seem to oc- 
cupy very large wall-like regions (Doroshkevich et al. 1996). The high density 
superclusters incorporate ~40% of all galaxies, and they have a typical diame- 
ter of ~30/i _1 Mpc (Hubble parameter h in units of 100 km s _1 Mpc -1 ). They 
surround huge underdense regions of diameter (50 - 70)/i~ 1 Mpc. 

The appearance of such a supercluster distribution is a typical outcome of 
the nonlinear gravitational instability as described by the Zel'dovich theory (e.g. 
Shandarin & Zel'dovich, 1989). The filamentary galaxy distribution has been 
characterized as forming a cosmic web with a sponge like structure (Bond et 
al. 1996). A quantitative analysis shows that the characteristics of very large 
structures and the evolution history are quite sensitive to the specific cosmo- 
logical model. In particular, the wall like structures are formed during a short 
evolution stage, and later they disintegrate into a system of high density clumps 
(Doroshkevich et al. 1997). The large scale matter distribution can be traced 
back to the initial peculiar gravitational potential (Madsen et al. 1998). Using 
statistical characteristics of the potential in different cosmological scenarios, the 
formation and evolution of the supercluster structures can be described quan- 
titatively (Demiahski & Doroshkevich 1998). These theoretical predictions and 
our numerical experiments can be used to constrain cosmological scenarios which 
give rise to the various peculiar potential fields. 
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We concentrate on a quantitative characterization of the geometrical prop- 
erties of the supercluster structures. We take up the attempts of Babul & Stark- 
man (1992) and Luo & Vishniac (1995) to quantify the linearity and planearity 
of the galaxy clustering with moments of the tensor of inertia. These statistics 
alone do not strongly discriminate between different cosmological models, since 
they characterize the general character of the anisotropic gravitational collapse. 
We perform a cluster analysis separately in high-density and low-density regions 
of the universe (the demarkation line lies about at an overdensity 10, see below), 
and we provide characteristics of the clusters and superclusters in dependence 
of the linking length used to define it. Our basic results lead to a remarkable 
discrimination between rich structure elements and low density regions in the 
universe. Recently, specific filament finders are investigated (Dave et al. 1997) 
which are able to follow the cosmic web. This statistic is strongly restricted to 
the small scale galaxy distribution. It discriminate marginally between different 
scenarios of structure formation, but no specific model is clearly favoured in the 
analysis of the CfA redshift survey. Here we employ the core sampling approach 
and the minimal spanning tree to define and to characterize the wall-like and the 
filamentary structures and to measure the parameters of their sizes, abundances 
and distributions. 

First we provide a description of the simulations for modeling the super- 
cluster distribution. Next we describe the cluster analysis and minimal spanning 
tree techniques used to analyse it. Then we employ the core sampling approach 
to characterize the supercluster structures in comparison to the LCRS. Finally 
we compare quantitative characteristics of the supercluster structures in depen- 
dence on theoretical predictions, and we discuss the results. More details can 
be found in Doroshkevich et al. (1998). 

2. Simulating the large scale structure 

We study the large scale matter distribution in N-body simulations using some 
of the promising cosmological scenarios presently under discussion. We use six 
cosmological models: A COBE normalized CDM model with density parameter 
£l m = 1 (SCDM) is taken as a comparison model despite of its well known 
difficulties. Alternatives based on CDM and a flat geometry of the universe 
include modifications of the primordial power spectrum by introducing a tilt 
oc /c ' 9 , or a break in the power spectrum at a certain scale, TCDM and BCDM, 
respectively. Both are inflation motivated (Gottlober et al. 1991,1994). Further 
we use open CDM models with fl m = 0.5 and VL m = 0.35, OCDM1 and OCDM2. 
Furthermore we take a model with f^o = 0.35 and a cosmological constant to 
provide spatial flatness, ACDM. All models are COBE normalized (Gorski et al. 
1994, Stompor 1995), the amplitude is characterized by the mass variance as on 
a scale of 8/i" 1 Mpc as given in Table 1. 

The simulations were run in boxes of very large comoving side length Lb ox = 
500/i -1 Mpc to include with high statistics the range of wave numbers k~ l f» (10 
- 30) /i _1 Mpc responsible for the supercluster formation. We use a PM code, 
described in more detail in Kates et al. (1995). We use N p = 300 3 particles 
in N ce u = 600 3 grid cells, which provides a force resolution ~ 0.9/i _1 Mpc, the 
mass resolution is given by m p in the Table 1. For later reference, we show also 
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Table 1. 



Parameters of simulations. 



model 



• m 



h 



10 n M 



km/s 
1375 
1293 
714 
550 
372 
913 



lo 

iMpc 



h- 1 



SCDM 
TCDM 
BCDM 
OCDM1 
OCDM2 

ACDM 



0.35 
0.35 



1 
1 
1 

0.5 



0.5 

0.5 

0.5 

0.6 

0.65 

0.7 



3.2 
3.2 
3.2 
1.6 
1.1 
1.1 



1.37 
1.25 
0.60 
0.74 
0.57 
1.12 



13.2 

13.2 

13.2 

22. 

29. 

26.9 



0.6 

0.57 

0.31 

0.14 

0.14 

0.2 



the velocity dispersion of the matter particles a ve i in the different simulations, 
which are scaled to the present time according to the linear perturbation theory. 
A length scale lo and a dimensionless constant tq characterize the typical scale 
of the supercluster structures, and their evolution stage as discussed below: 



Here k is the comoving wave number, the factor k stems from the input primor- 
dial Harrison-Zel'dovich spectrum, T(k) denotes the transfer function of pertur- 

1/3 

bations in the different scenarios, and k m i n = 2tt/Li )OX , and k max = k m i n N c ^ u 
gives the range of the spectrum realized in the simulations. 

In Retzlaff et al. (1998), we employed these simulations to discuss the power 
spectrum of the clustering of Abell clusters. There, model clusters were identi- 
fied with maxima in the smoothed density field. Now we use the dark matter 
particles to characterize the geometry of the superclustering. For comparison 
with the observed galaxy catalogues, we must allow for differences in the distri- 
bution of galaxies and dark matter. Here we start with the simple hypothesis 
that on large scales, the galaxy distribution follows the gross distribution of the 
dark matter. As shown by the mass variance of the different models, we must 
allow for a bias or an anti-bias in the small scale galaxy clustering to get the ob- 
served galaxy variance osfeal) = 1> an d to reproduce the observed galaxy- galaxy 
correlation function (Tucker et al. 1997). However, we show results only for the 
full dark matter distribution. They are independent of any additional assump- 
tions on prescriptions for the 'galaxy identification'. Effectively, we employ mass 
cuts in the cluster definitions which are equivalent to a biasing. For producing 
'galaxy catalogues' we also employed simple local bias prescriptions depending 
on the environmental density. Our conclusions are basically unchanged by such 
a biasing. The reason is that we study the geometrical properties of the matter 
distributions on very large scales which are independent on differences of the 
clustering of dark matter and galaxies in high density clumps. 

To illustrate the simulations, we show in Fig. 1 a cone diagram of a mock 
'galaxy sample' for the LCRS which is based on ACDM. We took the Schechter 
function of the LCRS-galaxies and assigned randomly luminosities to the dark 
matter mass particles according to the Schechter function characterizing the 
LCRS (Lin et al. 1996). Obviously, huge overdense regions are seen which 
match properties of the real observations. 
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3. Geometric characterization of large scale structures 

The cluster analysis is widely used in cosmology both to define virialised ob- 
jects, and for large linking lengths, to describe characteristics of the large scale 
matter distribution up to percolation, cf. e.g. Sahni & Coles (1995). It is 
closely related to the characterization of the galaxy distribution by the minimal 
spanning tree (MST, cf. Barrow et al. 1985). The MST represents an unique 
geometric construction in the galaxy distribution which is based solely on the 
galaxy positions. The simplest characterization of the MST is the number dis- 
tribution of the edge length, Wmst(1)- First, Barrow et al. (1985) constructed 
it for the CfA-catalogue. At small separations, this distribution characterizes 
the non-linear galaxy clustering, and on large scales, it provides a measure of 
the topology of the geometrical support of the galaxy distribution. For random 
particles situated on two-dimensional (walls) or one-dimensional structures (fil- 
aments) , we get logarithmic increments in the number of edge length of two and 
one, respectively. Parametrising the edge length distribution by 

Wmst(x) = -Wo ^ e-W, x = l/<l>, 

where I denotes the edge lengths of the MST and < I > its mean. Wq provides 
the normalization and the function Ff(x) is parametrised as 

F f {x) = (Pix Pl +(3 2 x P2 ) p - i . 

As mentioned, we are mainly interested in the effective power index p(x) = 
dlogFf/dlogx, characterizing the geometry of the galaxy distribution. 

Splitting off the MST at different linking lengths, we get standard friends- 
of- friends clusters with various linking lengths ru n k, the number of clusters 
N*t(runk) (the exponent t is used to indicate on the total number of clusters, we 
discuss separately also clusters with different multiplicity, singlets N^, doublets 
etc.): 

N'drunk) = / WmstWI. 



Table 2. Power law indices of cluster numbers and MST edge lengths 
in the ACDM model, TOT, RSE, and LDR denote total sample, rich 
structure elements, and low-dense regions, respectively; com in comov- 
ing and red in redshift space. 



sample 


fv 


< n p > 
h 3 Mpc~' 3 


P 3 


P* 


< Imst > 


pMST 










h-iMpc 




TOT-com 


1.0 


0.113 


0.71 ± 0.02 


0.59 ±0.03 


0.72 


0.60 ±0.02 


TOT-red 


1.0 


0.113 


0.76 ± 0.02 


0.85 ± 0.04 


0.78 


0.91 ± 0.03 


RSE-red 


0.45 


1.86 


1.82 ±0.02 


1.83 ± 0.03 


0.42 


1.60 ±0.04 


LDR-rcd 


0.55 


0.065 


1.01 ± 0.02 


1.11 ±0.02 


1.00 


1.11 ±0.03 



In Table 2 we provide the power law exponents of the edge length distri- 
butions of the MST (p MST ) and of the cluster number versus linking length (p s 
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for singlet clusters, p t for complete clusters) in simulated dark matter distribu- 
tions of the ACDM model as a typical example. We show the results in the 
physical space and, to allow a comparison with the data from observed redshift 
catalogues, in redshift space, were the galaxy positions are shifted along one 
axis according to their peculiar velocities. In addition we use a linking length 
kink of 1 /i _1 Mpc and a threshold multiplicity N thr = 200 for discriminating 
rich structure elements (RSE) with a minimum overdensity. Effectively this 
determines those particles in the simulations which mark 'superclusters' in the 
galaxy distribution. The remaining particles are separately analysed, they form 
low-density regions (LDR) which are characterized by field galaxies. The same 
discrimination was made by Doroshkevich et al. (1996) in the Las Campanas 
Redshift Survey. Similar as in the galaxy data, we identify a percentage / ~ 0.40 
in rich structure elements in some cosmological models. In the simulations we 
find power law exponents p s , p l , and p MST of about unity if analysed in physical 
space, i.e. the matter distribution is basically filamentary. In redshift space, the 
power laws for rich structure elements have exponents p l « 1.6 and p MST m 1.8, 
while again they are near unity in low-density regions. This difference is the 
main justification for the discrimination between different populations in the 
galaxy distribution. In redshift space the influence of the velocity dispersion 
erases the filamentary character of the matter distribution in RSE. Therefore 
we find an apparent particle distribution in huge wall- like structures. 

In the LCRS a power index p t ~ pmst ~ 1.7 has been found for RSE, and 
Pt ~ Pmst ~ 1 for LDR and for the total sample. This is comparable with our 
results in redshift space. A complex inner structure of RSE is also visually seen 
in the LCRS which resamples the galaxy distribution in the great wall of the 
CfA catalog (Ramella, Geller & Huchra 1992). 



4. Mean separation of filamentary and sheet-like components 

The core-sampling method of Buryak et al. (1994) allows us to discriminate the 
filamentary and sheet-like structure elements and to find quantitative charac- 
teristics of these structures, namely, the surface density of filaments, 07, that is 
the mean number of filaments crossing a randomly oriented unit area (in h~ 2 
Mpc 2 ), and the linear density of sheets, a s , that is the number of sheets crossing 
the unit length (in h^ 1 Mpc) of a random straight line. These parameters are 
equivalent to the mean separation between sheet-like structure elements, D s , 
and filaments, Dfi 

D s = l/a s , Df = aj 1/2 , 

i.e. these lengths represent the mean free path between sheet-like and filamen- 
tary structure elements. 

For our analysis, 196 cylinders with a radius of 1.7ft, _1 Mpc were prepared. 
The mean number of points within the cylindrical cores amounts ~(400 - 600). 
The analysis was performed for 16 values of the cylinder radius, 1.7/i _1 Mpc 
> f C yi > 0.7/i _1 Mpc. The mean separation of sheet-like elements, D s , and the 
surface density of the filamentary component cr/, are plotted in Fig. 2 as function 
of the fraction of matter which is randomly selected. We provide results only for 
the OCDM1 and ACDM models which are the only reasonable models as the next 
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Figure 2. Mean separation of the RSE, D s , and the surface density 
of filamentary component, <T/, vs. the matter fraction concentrated 
within the structures, /, in redshift space for the OCDM1 (dots) and 
ACDM (stars) models. 

section demonstrates. The range of values found in the LCRS (Doroshkevich et 
al. 1996) are marked by dotted lines. There is a clear signal from a wall-like 
supercluster population with a characteristic separation of Ds ~ (40 — 60)/i _1 
Mpc. But in contrast with results from the observations, D s increases slowly for 
small / < 0.6. This effect is probably caused by differences in the covering of the 
cores by the walls in the simulations, cf. also Ramella et al. (1992) and Buryak 
et al. (1994). On the other hand, the nearly constant D s (f) and a quick drop 
of cr/ with / proofs that a significant matter fraction (~ 0.4 — 0.5) is associated 
with the high dense sheet-like component, consistent with the observations. 

5. Quantitative description of the overdensity regions 

The discrimination of rich structure elements (RSE, or 'superclusters') and low- 
density regions (LDR populated by field galaxies) delivers no sharp criterion, 
we varied the linking length and the mass threshold for identifying structures. 
The results are shown in Fig. 3 and 4 for the various cosmological models. 
As discussed below, the discrimination is much more dependent on the type of 
cosmological model than on the exact boundary and criteria used. Therefore, 
the analysis delievers a test of cosmological models which is very robust. 

In Fig. 3 we show the dependence of the mean overdensity of clusters on the 
fraction of matter in the supercluster structures. The different symbols belong 
to different mass thresholds (here given in terms of the minimum number of 
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Figure 3. Overdensity Sodr m redshift space vs. matter fraction 
foDR for thresholds: N th = 100 (dots), N th = 200 (dashed line), N th = 
300 (stars), N t h = 500 (long-dashed line), N t h = 1000 (open squares). 
Dotted lines show the observed parameters from LCRS. 



particles per cluster). The dotted lines show the approximate characteristics of 
the LCRS (Doroshkevich, 1996). Obviously we get rich structures in the high 
density regions only if the model is well evolved, in particular, OCDM2 and 
BCDM cannot reproduce the degree of superclustering in the LCRS. 

In Fig. 4 we show the dependence of the cluster velocity dispersion along 
the minor axis in dependence on the fraction of matter within clusters. The 
velocity dispersion can be estimated from the data either using the extension 
of the superclusters along their smallest axis assuming virial equilibrium of the 
transverse galaxy motion. In this spirit, already Oort (1983) gave an estimate of 
the observed velocity dispersion a° bs ~ (350—400) km/s. The transverse velocity 
dispersion of superclusters in SCDM and TCDM is much to high, i.e. these 
models do not lead to the observed narrowness in the supercluster distribution. 
Also, ACDM shows a slight excess in the velocity dispersion along the smallest 
axis. The stablest superclustering is produced in OCDMI. It is remarkable that 
the OCDM2 models has a too low velocity dispersion and a too low concentration 
of matter in superclusters. This model is simply underevolved. 

For SCDM and ACDM we show a few evolution steps. It is obvious that 
the supercluster structures are built up quite recently in both models. Also 
the transverse velocity dispersion develops strongly due to the nonlinear grav- 
itational clustering and the disintegration of the wall-like structures in huge 
filaments and high-density clumps with high and isotropic velocity dispersion. 
The strong evolution of the superclustering up to redshift z = 3 has to be 
compared with some observational evidence of large matter inhomogeneities at 
redshifts z « 2 — 3 (Connolly et al 1996, Cristiani et al 1996, Williger et al. 
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Figure 4. Velocity dispersion a vt vs. matter fraction concentrated 
within rich structure elements Fuse for the same richness thresholds 
as in Fig. 3. Dotted lines show the observed parameters in the LCRS. 



1996). Compared to our simulation results, such structures must be very rare at 
high redshifts, or there must be a strongly evolving bias of the galaxy formation 
up to these redshifts. 

We use the eigenvalues of the tensor of inertia to estimate the degree of 
anisotropic gravitational collapse in forming the supercluster structures. The 
mean length of the rich structure elements amounts L\ ~ 16, 17, and 24/i _1 Mpc 
for SCDM, OCDM1 and ACDM, respectively. Along the other axes we find a 
collapse by a factor of two and (4 — 5), respectively. This is consistent with the 
mean overdensity of 10 for the RSE, and we expect an initial extension of the 
superclusters in the simulations of Lf^ t « (20 — 25)/i _1 Mpc (comoving coordi- 
nates), i.e. no collapse of the structures along the longest axis. This estimate 
enters in the following theoretical estimates. 



6. Theoretical estimates of supercluster parameters 

The reproduction of the observed characteristics of the RSE in simulations with 
a standard CDM-like power spectrum verifies that the observed structure may 
have formed during the nonlinear evolution of small initial perturbations, and 
that its characteristics can be derived from the parameters of the initial power 
spectrum of Gaussian fluctuations, and a specific cosmological model. Demiahski 
& Doroshkevich (1998) delivered an approximate theoretical description based 
on Zel'dovichs nonlinear theory of gravitational instability. The comparison of 
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the initial size of structures, L init , the wall separations, D s , and the mean ve- 
locity dispersions in rich supercluster structures, a u , and the velocity dispersion 
along the smalles axis of the superclusters, 0-3, are most interesting. They are 
connected to the length scale Iq and evolution parameter tq of the power spec- 
trum: 

hf lt ~ 2.9r Zo, 
and for the SCDM, OCDM1 and ACDM models we have 

Lf nit m 23fc -1 Mpc, w 9fc -1 Mpc, « 16h -1 Mpc, 

respectively. For the OCDM1 and ACDM models Lf nit w (0.6 - 0.7)LJ*g is 
found. The difference to our numerical results is caused partly by the mass 
threshold used, and partly by the estimates of the cluster boundary. In the 
SCDM model we get Lf it ps 1.5LJ£™, it is strongly influenced by the small scale 
clustering. 

A theoretical estimate to the mean wall separation for the mean rich super- 
clusters is given by 

Df w 6 • IqV^, 
and for the three models above we have 

Df w 56/i _1 Mpc, w 49/T 1 Mpc, w 72/T 1 Mpc, 

respectively. These values are consistent with the values of D s found with the 
core-sampling method. 

The theoretical model shows that for Gaussian perturbations the distribu- 
tion of the random velocity of structure elements is also Gaussian with a velocity 
dispersion 

af « y/3H (/3 - l)kr = (P - l)a pvd 
where the parameter (3 characterizes approximately the evolutionary stage of 
perturbations, for SCDM = 2, otherwise 

(l + 4J2 m )/(l + 1.50 m ) ps 1.7 for OCDM1 
/3 » (1 + 3.4^ m )/(l + l.2n m ) w 1.54 for ACDM. 

For rich structure elements, af 1 is lower by a factor ~ ^JJqJl^^ = (2.9ro) -1 ^ 2 . 
Then for SCDM, OCDM1, and ACDM, we have 

a u » 950 km/s , 385 km/s , 490 km/s , 

respectively. These estimates are consistent to the simulations. 

A theoretical estimate for the inner velocity dispersion in RSE along the 
smallest axis 03 takes into account only the matter infall into RSE, and it ignores 
the disruption of the pancakes into high-density clumps. We have 

and for a typical primordial size, L*raif gi ven above, 

af ~ 0.5(/3 - l)a pve i ^ 0.5a u . 

This expression should be taken as a lower limit for the inner velocity dispersion, 
and, really, for all models under consideration, we have af ~ (0.5 - 0.6) a 3 , i.e. 
the influence of small scale clustering on the properties of RSE is very important. 



10 



7. Discussion 



In our analysis the properties of simulated spatial matter distributions were 
studied for six cosmological models with CDM-like power spectra. The simula- 
tions were performed in large boxes in order to reproduce correctly the mutual 
interaction of large and small scale perturbations, and to obtain a representative 
sample of wall- like RSE. The broad set of considered cosmological models allows 
us to reveal the influence of main cosmological parameters on the formation and 
evolution of the wall-like RSE, and to discriminate between these models. Our 
results show that the used methods are effective, and they yield a description 
of the spatial matter distribution as a whole and, in particular, the character- 
istics of the supercluster distribution. The main results of our analysis can be 
summarized as following: 

• The phenomenon of the strong matter concentration of galaxies in wall- 
like rich structures can be reproduced with standard COBE normalized 
CDM-like power spectra for suitable cosmological models and simulation 
parameters. An essential fraction of DM, f rse ~ 0.4, is compressed non- 
linearly on the scales ~ (15 - 20) /i _1 Mpc that is less then the mean 
separation of these RSE by a factor of ~ 2 - 3. 

• The comparison of observed and simulated parameters of the wall-like 
supercluster distribution allows us to discriminate between different cos- 
mological models. Only the low density models, the ACDM model with 
tt m h « 0.15 - 0.25, or the OCDM1 model with Q m h « 0.25 - 0.35 can re- 
produce the observed concentration of galaxies in rich superclusters. Prob- 
ably, similar results can be obtained for a model with a mixture of cold 
and hot dark matter. A large scale bias between the spatial distribution 
of DM and galaxies can also improve the predictions of some models. 

• Our results underline theoretical expectations concerning the epoch of the 
supercluster formation. At z = 1 the fraction of matter accumulated by 
superclusters with the chosen richness and overdensity drops by a factor 
~ 2, and at z = 3, it becomes very small. 

• The simulations reproduces both the wall-like supercluster distribution and 
the filamentary component of the low density regions. Each component 
represents an essential fraction of matter in the universe, and it is equally 
important for the description of the joint network structure in the spatial 
matter distribution. 

Taken together only the low density models, best OCDM1, but also ACDM, 
can reproduce the large sheetlike matter distribution observed in recent redshift 
surveys. 
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